Landau Expansion for the Critical Point of a Polydisperse System 
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The effect of polydispersity on the phase diagram of a 
simple binary mixture is to split the binodal curve into cloud 
and shadow curves that cross at the critical point (which, in 
general, is not at the maximum of either curve). Recent the- 
ories of polydispersity have shown, in favorable cases, how to 
project the (infinite-dimensional) free energy of the polydis- 
perse system into a low dimensional space of 'moment densi- 
ties'. We address here the issue of how to construct a Landau 
expansion from the projected free energy. For the simplest 
case where the excess free energy depends on one moment 
density pi (this includes Flory Huggins theory for length- 
polydisperse chains) we show that the minimal expansion re- 
mains quartic in pi but nonetheless has seven independent 
coefficients, not two. When this expansion is handled cor- 
rectly all the peculiar qualitative features of the polydisperse 
phase diagram are recovered, as are the slopes of the cloud 
and shadow curves, and the curvature of the cloud. However, 
to correctly predict the curvature of the shadow, certain fifth 
order terms must be retained. We also consider the phase 
diagram on the temperature-pressure plane, where the co- 
existence line broadens into a region. In general, the critical 
point lies between the point of maximum temperature and the 
point of maximum pressure on the perimeter of this region. 
This behavior is again captured by the Landau expansion. 

PACS numbers: 



I. INTRODUCTION 

Many systems, such as colloidal suspensions, and poly- 
mer blends and solutions, contain particles with an effec- 
tively continuous distributions of an attribute (cr, say) 
such as size, charge, or chemical character. Such sys- 
tems are termed polydisperse [1]. Their phase equilibria 
can become complicated, but may exhibit relatively sim- 
ple features including coexistence of two liquid phases 
of different composition. In some cases the tendency to 
phase separate is inherited from the underlying monodis- 
perse physics (e.g. attractive colloids); in others it stems 
only from polydispersity (e.g. fractionation among chains 
with a spectrum of chemical compositions). 

In a polydisperse system undergoing fluid-fluid demix- 
ing, the phase diagram looks rather different from that of 
a simple binary mixture even if the underlying physics is 
similar. It is common experimentally to draw the phase 
diagram on a composition/temperature plane; one effect 
of polydispersity is then to split the binodal curve into 



cloud and shadow curves that cross at the critical point. 
(The latter is not at the maximum of either curve.) This 
is because, when a system is cooled to the point where it 
becomes globally unstable to phase separation (the cloud 
point), the incipient new phase (its shadow) differs from 
the old one not only in the number of particles per unit 
volume, as usual, but in its cr-distribution. Because of 
this difference, the incipient new phase does not itself 
represent a point on the cloud curve. (However, it does 
lie on the cloud curve of a different system whose over- 
all size distribution coincides with that of the incipient, 
shadow phase). 

An interesting question is whether all this behavior 
can be captured within a Landau expansion. There are 
at least two reasons to consider this. Firstly, concep- 
tual (as well as practical) progress has been made re- 
cently in understanding phase equilibria in polydisperse 
systems within a certain class of mean-field models: those 
in which the excess free energy depends on the a distribu- 
tion only through a finite set of moments [1] . A promising 
method involves projecting the free energy onto a rela- 
tively small number of densities; when these are chosen 
correctly, the cloud, shadow, and spinodal curves, and 
with them any critical points, are found without approx- 
imation beyond those used to set up the mean-field free 
energy in the first place. 

Due to their mean-field character and restricted set of 
relevant densities, such projection techniques should be 
reconcilable with a Landau expansion near the critical 
point - yet the unusual topology of the phase diagram 
appears, at first sight, to preclude this. In particular, 
the usual adoption of a quartic expansion in one den- 
sity (with coefficients dependent only on temperature) is 
bound to exhibit a critical point at the apex of a bin- 
odal curve. We show below that a more careful Landau 
expansion does not suffer from this restriction, and in- 
stead reproduces the correct, more complicated behavior 
outlined above. 

A second reason to study the Landau expansion is in 
the hope of going beyond mean field theory. Experi- 
mental values for the critical exponents of a polydisperse 
system [2] deviate from those of the monodisperse case 
and are consistent with Fisher-renormalised exponents 
[3]. However, to our knowledge, there is no explicit renor- 
malization group calculation for a polydisperse system. 
Also, quite apart from calculating critical indices, it is 
important for many purposes (see e.g. [4]) to properly 
classify and identify the scaling fields near the critical 
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point, and this task cannot really be begun without a 
consistent description at the Landau level. Hence the 

Landau expansion presented below represents a first step 
towards a fuller understanding of polydisperse critical be- 
havior beyond mean field level. 



II. MONODISPERSE SYSTEM 

In this section we review briefly the classical Landau 
expansion of the free energy around the critical point 
for a simple monodisperse system. This will facilitate a 
comparison with the polydisperse case and establish the 
notation. 

In a canonical ensemble, the free energy T is defined 
as a function of the temperature T, the number of parti- 
cles N and the volume V . In the thermodynamic limit, 
however, we can define a reduced free energy (density) as 



HN,T,V) 



(1) 



which is a function only of the density of particles p = 
N/V and the temperature. At fixed temperature, any 
terms that arc constant or linear in density are irrelevant 
for obtaining the phase equilibrium of the system (which 
is our main aim) since they can be assimilated into shifts 
of the entropy and of the chemical potential. 



A. Critical Point 

Within the mean field approximation, which is where 
the Landau expansion holds, the critical point must obey 
one condition common to all points on the spinodal curve. 



= 0, 



(2) 



and a second condition particular to the critical point 



9p3 



0, 



(3) 



where |c means evaluated at the critical point {p^jTc). 



B. Landau Expansion 

Once the critical point is determined, the expansion 
follows as usual: 
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(4) 



where 6p = p—p°. We forget about the constant and lin- 
ear terms since, as described above, they can be viewed as 
(temperature dependent) shifts in the entropy and chem- 
ical potential with no effect on phase behavior. 

The coefficients of the expansion depend on the tem- 
perature and can be expanded themselves about the crit- 
ical temperature. Using Eqs. (2) and (3), and keeping the 
lowest order term for each coefficient, we obtain: 



1 d^F 



STSp^ 



1 d^F 



3! dp^dT 



3^1 d^F 
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Sp' + .... 



(5) 



where ST = T - T^. Since 5T ~ 5p^ (see Eq. (7) below) 
the term in 5T5p^ is of higher order and can be neglected. 
Therefore, the lowest order expansion of the free energy 
has the simple form: 



AdTSp^ + C Sp'^, 



(6) 



for appropriate A and C. The coexisting phases obey the 
standard result 



(7) 



giving the well-known parabolic shape of the phase dia- 
gram around the critical point for monodisperse systems 
(always within the mean field approximation). 



III. POLYDISPERSE SYSTEM 

In the polydisperse case, the particles are not all identi- 
cal but differ in certain scalar attributes (particle radius, 
charge, chain length...) which can be tagged with indices 
ai, (72, . . .. For simplicity, we focus on systems with just 
one polydisperse index a, and write the corresponding 
density of particles as p(cr). In this case, the free energy 
can be written as two parts, an ideal term and an excess 
term, as follows [5]: 



F[p]= jda p(a){log(p(a))-l} 



(8) 



While the ideal part only depends on p{a) as shown, 
the excess part Fex is, in general, an arbitrary functional 
of p{a) and a function of T. Here, in order to explore 
how polydispersity affects the Landau expansion around 
the critical point in a tractable manner, we restrict our- 
selves to cases where the functional dependence of the 
excess part of the free energy on the density distribution 
p{a) is through a single 'generalized moment' [6] of the 
distribution: 



F[p]= jda p(a){log(p(a))-l} + /(pi,T), 



(9) 



where 
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pi= Jda p{a)wi{a) (10) 

with wi{a) an arbitrary weight function. Although this 
assumption seems restrictive, it is true; for some im- 
portant models including the Flory-Huggins theory of 
length-polydisperse polymer solutions and blends [6]. 
The models for which the excess part of the free energy 
depends on the density p{a) through a finite set of mo- 
ments of this distribution, for which Eq. (9) is a particular 
example, are usually termed truncatable models. 

A. The Parent Distribution Function 

Polydispersity typically appears in the process of 
chemical manufacture of the particles of the system, so 
that the relative proportion of particles with different a 
is fixed thereafter. The actual density of each species 
then depends on how much solvent is present. We define 
a parent distribution P{a) as 

P{a) = Po (11) 

where n{a) is a normalized distribution (fixed by the ini- 
tial preparation of the system) and Pq is the overall den- 
sity of particles, which varies as one moves along a 'dilu- 
tion line' through the phase diagram. All other moments 
of the parent are fixed by Pq; for example. 

Pi = jda P{a)wi{a) = Pq ni (12) 

where rii is fixed at synthesis via the parent: 

rii = Jda n{a)wi{a). (13) 

Although the parent distribution is globally fixed, the 
system can separate into coexisting phases with a distri- 
bution p{a) that is locally different from P{a). If there 
are two such phases a and /3, with respective distribu- 
tions p"(cr) and p^{cr), the parent distribution is recov- 
ered as the volumetric average: 

x^p'^ia) + x^p^ia) = Pia), (14) 

where Vx" and Vx^ are the volumes of phases a and 
(3, respectively. Note however that at a phase boundary 
which marks the onset of phase separation, the volume 
fraction of the incipient (shadow) phase is negligible and, 
therefore, the cr-distribution of the original (cloud) phase 
must coincide with the parent distribution. 

B. Projection Method 

Although for a given n{a) we require only one density 
(such as pi), alongside temperature, to specify a point on 



the phase diagram, the corresponding free energy must 
depend on the full density distribution p{a) = Pon{a). 
Accordingly it is not correct to substitute the parent dis- 
tribution in a functional like Eq. (9) to obtain a one- 
variable free energy F{Po) (or equivalently, F{pi)). In 
fact such a substitution would describe a case where 
the normalized distribution n{a) were constrained to be 
equal in all coexisting phases - which is not the physical 
situation. Instead, in the case of Eq. (9), the equality 
of chemical potentials dictates that the a-distribution of 
any cloud phase (which coexists, by definition, with an 
unperturbed parent) must obey 

p{a) = P{a) e^'"'^''\ (15) 

for an appropriate value of A. This distribution is not 
proportional to the parent unless A = (when the phases 
are identical). 

However, under certain conditions there is a procedure 
to transform the free energy functional into an effective 
free energy function. The conditions are first that (as 
in Eq. (9)) the excess free energy / only depends on the 
distribution p{a) through a finite set of moments (defined 
as in Eq. (10)); and second that we are only interested 
in phase boundaries (the cloud and the shadow curves) 
and the spinodals. The procedure, called the projection 
method [6,1], involves several steps briefly described here, 
and leads to a free energy function dependent only on the 
moments required to specify the excess free energy. 

First, considers the free energy 

F[p\= ydap(a)|log(|^)-l|+/(pi,T), (16) 

which, for our purposes, is equivalent to Eq. (9), since 
the inclusion of the parent as a denominator in the loga- 
rithm represents a linear term (which can be considered 
as a mere shift of the chemical potentials). But note that, 
for that equivalence to hold, the parent -P(c), including 
Po) must be kept fixed as the density distribution p{a) 
varies: this observation is crucial to the Landau expan- 
sion developed below. 

The projection method follows by minimising the first 
(ideal) term of this free energy subject to the constraint 
of fixed values for the moment (s) upon which the excess 
free energy depends. (The only such moment is pi in our 
case.) This yields 

p{a)=P{a)e^''"^^''\ (17) 

where Ai is the Lagrange multiplier used to constrain the 
minimization, determined by p\ through 

pi = jda P{a) wi (a) e ""^ . (18) 

Since both the parent phase (the cloud) and the coexist- 
ing phase (the shadow, Eq. (15)) obey Eq. (17) for differ- 
ent values of Ai, the method contains no approximation 
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when it comes to finding the phase boundaries. However, 
if the cr-distribution of any of the phases adopts a form 

different from (17), the method produces only approxi- 
mate results. This occurs whenever the system separates 
into finite volumes of more than one phase, that is, every- 
where within the cloud curve. (In Eq. (14), the projection 
method is inexact when Xa and are both nonzero.) 
The free energy obtained by this method is [6] 



F{pi,T) = Ai pi - po + f{pi,T), 



where 



(19) 



(20) 



This free energy depends on only one density variable, 
since both Ai and po are given by pi. The coexisting 
phases are obtained from this free energy function via the 
common tangent method, as in the monodisperse case. 



C. Dilution Line 

As mentioned above, the projected free energy re- 
quires knowledge of the parent distribution function 
P((7), whose normalization Pq varies as one moves along 
a dilution line on the phase diagram. For this reason, 
the projected free energy F{pi,T) depends parametri- 
cally on Pg, through Ai(pi,Po) and po(Ai,Po); see Eqs. 
(18) and (20), respectively. This dependence is paramet- 
ric because, in the process of calculating the coexistence 
between phases, Pq must remain fixed, so as to maintain 
the equivalence of Eq. (9) and (16). 

Suppose that, for fixed Pq and T, we obtain by this 
method a coexistence between phases a and /?, with mo- 
ment densities p" and pf , respectively. This determines 
the value of A", A^, Pq, and pg through Eqs. (18) and 
(20) which, in turn, determines the density distribution of 
the two phases p"(cr) and p^{a) via Eq. (17). In general, 
neither of these distributions is the parent distribution. 

This is consistent with the fact that the parent distri- 
bution does not coexist with a second phase for arbitrary 
combinations of Pq and T; this happens only on the cloud 
curve. To find the cloud curve, we combine the coexis- 
tence conditions derived above with the equation of the 
dilution line, so that one of the coexisting phases {a say) 
is the parent distribution; this implies that p"{a)—P{a), 
A« = 0, p^ = Po, andp? = Poni (or pf=p^ni). 



D. Critical Point 

The critical point of the polydisperse system lies at the 
intersection of the cloud and shadow curves and is, there- 
fore, correctly described (given Eq. (9)) by the projected 
free energy. The equations to determine it are analogous 
to those of the monodisperse case, 



and 



dpi 



g3p 



dpi 



= 



= 0, 



(21) 



(22) 



where \c means evaluated at the critical point 
{pIjTc, Pq). Note that these two equations have three 
unknowns, namely pj, T^, and Pq, and must be comple- 
mented by the dilution line constraint, pf =Po^ni. 



IV. LANDAU EXPANSION 

The projected free energy P(pi, T; Pq) can now be ex- 
panded about the critical point. 



P« - 
2 



1 d^F 



dpi 



Sp'' + 
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3! 



dpi 
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4! 



d-^F 



dpi 



Sp' 



(23) 



Pi,T;Po 



where 5p = pi — p^. As before, the constant and linear 
terms have been dismissed. 

Now we must expand the coefficients of this expansion 
around the critical values not only of T, as we did in the 
monodisperse case, but also of Pq. One might think that, 
after doing this, the free energy could be written, to the 
lowest order, as 



{A1ST + A2SP0) 6p'' + CSp\ 



(24) 



where 5Pq = Pq — Pq. If so, the term in brackets should 
be of order 6p^ to balance the term of order (5p^. How- 
ever, the second term in the brackets is actually linear in 
Sp, as dictated by the dilution line SPq = 5p/ni. These 
two facts can only be obeyed if the first term in the brack- 
ets behaves like A^ST « —A25p/ni + O {Sp'^) (instead 
of ST ^ Sp^, as in the monodisperse case). This corre- 
sponds to the fact that, generically, the critical point is 
no longer at the top of the coexistence curve. Therefore, 
to maintain consistently the order of the approximation, 
the correct lowest order expansion of the free energy must 
contain all terms of order ^p*: 



dpldT 
1 d^F 



6T + 



d^F 



dpldPo 



2 dpldP^ 



<5Po^ + 



<5Po + 

d^F 



1 a^F 



2 dpldT'^ 



dpi dTdPo 



6T SPq 



1 
3! 



d^F 



dpidT 



ST + 



d^F 



dpfdPo 
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SPo 



ST^+ 



Sp^ + 



Sp^ + 



4! dpi 

To continue, let us define the functions 



<5/ 



(25) 
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Pm = jda P{a) {wi{a 



(26) 



for any m = 0,1,2,..., which are an expanded set of gen- 
eralised moments consistent with the previous definitions 
Eq. (18) and (20). These functions depend on pi and Pq 
and, indeed, they obey: 



and 



dPo 



dpm _ Pm+1 

dpi P2 



1 / Pi Pm+1 

[Pm 

Po V P2 



(27) 



(28) 



Very importantly, these equations imply that the 
derivatives of any (generalised) moment with respect to 
pi gives rise to functions of moments (Eq. (27)), and the 
derivatives of any moment with respect to Pq must vanish 
in the monodisperse limit (as does Eq. (28)). 

With these definitions, we can write the second deriva- 
tive of the free energy, from which all the coefficients of 
the expansion (25) can be obtained, as: 



j_^aV(pi,T) 

P2 



dpi^ 



(29) 



Now note that, since T and Pq appear in disjunct parts of 
the free energy (as in the original free energy, Eq. (19)), 
all crossed terms must vanish: 



dTdPo 



= 0. 



(30) 



Accordingly, the coefficients of the expansion can be di- 
vided into three different groups: those including only 
derivatives with respect to pi, those including deriva- 
tives with respect to T, and those including derivatives 
with respect to Pq. While the first group of coeSicients 
originates from the full free energy, the second originates 
from the excess part of the free energy ('the interac- 
tions'), and the third group arises from the ideal part 
of the free energy ('the entropy'). These three groups of 
coefficients have a distinct dependence on polydispersity. 
In particular, all the coefficients stemming from the ideal 
part (those multiplying any power of 6Po) must vanish 
in the monodisperse limit [7], hence recovering the ordi- 
nary Landau expansion, Eq. (6) (plus some terms e.g. in 
5T^6p^ that are higher order corrections in this limit). 

In brief, the correct lowest order expansion of the free 
energy, retaining all terms to fourth order in density, is 
therefore 

F ^ [Ai 5T + A2 SPo + Aa ST^ + A4 5Pa^) 5p^ + 

{BidT + B2 6Po) dp'^ + CoSp^ , (31) 

where the coefficients of terms involving 5T (those with 
odd indices: Ai, As and Bi) only depend on derivatives 
of the excess free energy /, and the coefficients of terms 



involving SPq {A2, A4 and B2) stem solely from the ideal 
entropy term of the system and not the interactions. The 
latter vanish in the monodisperse limit; for example. 



1 



2P/^ 



n2^ 



1 

n2 



(32) 



where ni, n2 ■ ■ ■ arc generalized moments of the (normal- 
ized) parent distribution. 



Jda n{a) (wi(cr))^ 



(33) 



which do not depend on any thermodynamic variable. 

It is interesting to point out here that the ultimate rea- 
son why the critical point is not on the top of the cloud 
curve is that the squared term of the free energy expan- 
sion at the critical point behaves effectively as ST^ Sp^ 
instead of the usual 6T6p^. 

We stress that although the coefficients A2, A4 and 
B2 arc small for narrow polydispcrsities, the expansion 
(31) is asymptotically exact around the critical point for 
arbitrary (not necessarily narrow) polydisperse distribu- 
tions. As ever, 'exact' means exact within a mean-field 
free energy of the form Eq. (9). 



A. Cloud Curve 

In order to obtain an expression for the cloud curve, 
we can write the free energy expansion (31) in terms of 
a shifted density 



Sr] = 6p + 



Bi ()T + B2 i)Py) 



(34) 



obtaining an asymptotic expression of the form 

F^ASif + Co St]'^, (35) 

where 



A = Ai6T + A2 6Pn + A;, ST^ + A4 5Po 
3{B,dT + B2dPi>f 



Co 



(36) 



Starting from this simple expression, the coexisting 
states follow immediately as 



^77 w ± 



A 



(37) 



These two coexisting densities correspond to distribution 
functions given by Eq. (17), neither of which is in general 
the parent distribution function P{(t). In order to obtain 
the appropriate cloud curve, we must therefore impose 
the dilution line constraint, 6p'^'^ — ni SPq, which leads 
to the following expression for the cloud curve: 



STcL « aoL Sp + pcL Sp" + C>(^p^), 



(38) 
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where the slope of the cloud curve at the critical point is 
given by 



and its curvature by 

Co 



(1- A) -3 



Aa 



AlA, 



Ami Alnr 



where 



A = 



2 Ai Co ni 



(39) 



(40) 



(41) 



As predicted (preceding Eq. (25) above), we have thus 
established that Ai 6T + A2 5Po ~ 6p^, with exact can- 
celation of the terms linear in 5 p. This shows the consis- 
tency of the approximation and proves that none of the 
terms in Eq. (25) can be discarded in the expansion. 

One consequence of these results concerns the geomet- 
ric location of the critical point along the cloud curve. For 
the critical point to remain at the apex of the cloud curve 
(which is where it resides in the monodisperse limit), the 
slope of the cloud curve at the critical point act: Eq- 
(39), must vanish; this requires A2 = 0. According to 
Eq. (32), this requires very particular combinations of 
the parent distribution and weight function wi {a) . If the 
latter only adopts positive values, A2 is strictly positive 
except in the monodisperse limit so that the critical point 
will move away from the top of the cloud curve even for 
infinitesimal polydispersity [8] . Since the argument orig- 
inates from the entropic part of the free energy, it may 
hold more generally than for systems obeying Eq. (9) 
[9,10]. 

For monodisperse systems, the slope act vanishes, 
while tends to —2Co/Ai, thus recovering the classi- 
cal scenario that leads to Eq. (7): the cloud curve forms 
the binodal curve with a quadratic maximum where the 
critical point resides. 



B. Shadow Curve 

The shadow curve does not, of course, obey the dilu- 
tion line constraint, and therefore is not obtainable direct 
from Eq. (37). Nonetheless we expect 



STsH ~ asH Sp + PsH Sp'' + 0{Sp^). 



(42) 



To find asH, we may use Eq. (37) to write a relation 
between the cloud and the shadow curves. 



(5?7°^ + St]^" « 0. 
This yields, upon substitution of Eq. (34), 

acL 



oisH = - 



1 + A 



(43) 



(44) 



Interestingly though, Eq. (43) does not provide the 
correct curvature of the shadow curve, ^sh, as can easily 
be checked numerically for the example considered below 
(Section V). This is because the curvature of the shadow 
in fact depends on some higher-order terms (of order dp^) 
that were not included in Eq. (31). The fifth-order terms 
are 



S3 ST-" 
4 



BiSPo^) dp^ 



AF = {A5 6T^ + Aq 6Po^) Sp^ 

+ {Ci ST + C2 SPo) Sp" + Do Sp'' , (45) 

Including these (fifth-order) terms in the expansion, Eq. 
(31), transforms Eq. (43) into 



Sv""- + Sr,'" « 
(Ai ST + A2 5Po + Aj ST^ + ^4 SPq'') Do 
2 (Co + Ci5T + C2 SPof 

from which we finally obtain the expresion 



(46) 



+ 0{Srf), 



( 1 + A )'^ 



/3c 



1 + 



1 



2C>i 



(B2C0 - A2D0) 



+ 



2Clnl 



{B2C2 - BiCo + AiDo) 



+ 



(B1C2 + B2C1) 



OL, 



+:7^iBiC^ - B3C0 + A3D0) 



2 CI 

In the monodisperse limit, analogously to the cloud 
curve, the slope agn vanishes and pgn tends to —2Co/Ai, 
thus forming the binodal curve by the merger of cloud 
and shadow. 



C. Consequences of the Expansion 

Note the progression in complexity of the coefficients 
needed in the expansion to compute the slopes, and then 
curvatures, of the cloud and shadow at the critical point. 
The slope of the cloud curve act depends on just two 
coefficients of the expansion {Ai and A2); the slope of the 
shadow involves three extra coefficients {Bi, B2 and Co); 
the curvature of the cloud at the critical point requires 
knowledge of all seven coefficients of the quartic Landau 
expansion, Eq. (31); that of the shadow also involves a 
selection of the fifth order terms (but not all of them): 
Table 1 shows which terms in the Landau expansion must 
be retained to correctly predict each of these parameters. 

Although these results are obtained for a one-moment 
free-energy, Eq. (9), they involve higher moments of the 
parent distribution [6], as detailed in Table 1. Note also 
that the critical point and the spinodals depend only on 
a reduced set of moments; this holds to all orders in the 
Landau expansion. In contrast, order by order, the cloud 
and shadow curves involve progressively all the moments 
of the parent distribution, through Eqs. (27) and (28). 



(47) 
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The results given above were obtained by the projec- 
tion method, but since this gives the exact cloud and 

shadow curves, any other method for finding Ucl, Pcl ■ ■ ■ 
within mean field theory must produce, necessarily, ex- 
pressions equivalent to Eqs. (39,40,44,47). 



D. Phase Diagram in Temperature-Pressure Plane 

In the monodispcrsc case, the phase equilibria show up 
in the T—p phase diagram as a line of coexistence, ending 
at the critical point. If polydispersity is included, this line 
now broadens and adopts a banana-like shape of finite 
area, containing the values of temperature and pressure 
for which the system undergoes phase separation. In this 
case, the critical point appears on the perimeter of the 
coexistence region. 

We now can calculate the shape of the coexistence re- 
gion to the same order of approximation as we did with 
the cloud and the shadow. After projection of the free 
energy to give Eq. (19), the pressure of the system can 
be written in the familiar form [6] 



P 



dF 



Pi - F, 



ksT dpi 
which allows us to expand 6p = p — Pc as 



Sp 



ai ST + a2 SPo + 6i 5T^ + 62 5Po^ 



(48) 



(49) 



where (just as occurred in Eq. (31)) the coefficients of 
terms involving ST (those with odd indices: ai and 61) 
only depend on derivatives of the excess free energy /, 
and the coefficients of terms involving SPq (02 and 62) 
derive solely from the ideal entropy of the system (not 
the interactions) and vanish in the monodisperse limit. 
For instance, 



a2 



n 



^-1 (<0). 
n2 



(50) 



Substitution of the approximation derived in Section 
IV A for the cloud curve yields 



Up Sp + 0p Sp^, 



where 



and 



Qp = oi acL H 

m 



Pp = ai PcL + bi oici?' -h 



2 • 



(51) 



(52) 



(53) 



A generic consequence of these results concerns the lo- 
cation of the critical point along the perimeter of the 
coexistence region on the T — P phase diagram. It is 



clear on geometrical grounds that the critical point can- 
not be located at the point of maximum temperature 
and maximum pressure simultaneously, unless the coex- 
istence region shows a kink at that point. Such a kink 
is not predicted by the Landau expansion, which gives 
smooth behavior everywhere; nor have we seen evidence 
for it in any of the models (or experiments) studied to 
date. Assuming these two turning points are indeed dis- 
tinct, one can ask whether the critical point generically 
lies at one, at the other, or neither. We argued already 
above that generically, for systems obeying Eq. (9), the 
critical point is not at the maximum of T since this would 
require Ai to vanish (Eq. 32). Likewise, for it to lie at the 
pressure maximum would require ttp = in Eq. (52), for 
which a fine tuning of the parameters appears necessary 
[11,10]. We conclude that generically it lies at neither of 
these points; indeed, as shown in the following example, 
it lies between them. 



V. EXAMPLE 

In this Section, we illustrate the results of the Landau 
expansion for the cloud and shadow curves of a particular 
system: an incompressible blend of two homopolymers. 
One of the polymers is considered polydisperse in length, 
with number densities p((j) where u is the number of 
monomers on a chain (effectively a continuous parame- 
ter). The other polymer (or solvent) is monodisperse and 
has A''' monomers. For this system, the free energy can 
be easily obtained from the Flory-Huggins theory and re- 
covers the form (9). Specifically, the excess free energy 
reads [6] 



/(pi,x) = ^-^log 



N' 



+ XPi(l-/Ci), (54) 



where x is the Flory parameter, which plays the role 
of an inverse temperature. The weight function for the 
moment p\ is Wi(cr) = a (in the usual polymer jargon, 
p\ is the monomer volume fraction ^). 

As a normalised parent, we consider in this example a 
uniform distribution 



1 



2V3AAf 

otherwise 



if \a -N\< V3AN 



(55) 



where N =ni— {a) is the mean polymer size, and AN = 
((cr^) — (ct)^)^/^ is the standard deviation. 

We plot in Figures 1 and 2 the resulting behavior for 
the case where the polymer solvent has N' = 10 monomers 
and the distribution of sizes of the solute follows Eq. (55) 
with N = 1000 and AN/N = 0.2. (This corresponds, 
in polymer jargon, to a polydispersity index Myj/Mn = 
{a^)/{ar = 1.4.) 

We see from these plots that the Landau expansion, 
Eq. (31), captures correctly all the essential features of 
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the full mean-field model near the critical region, just as 
it does for monodisperse systems. In this example, the 
Landau predictions cross the cloud and shadow curves at 
the critical point (although the slopes and curvatures are 
correct). This behavior is inherited from the monodis- 
perse case, and is presumably connected with an error in 
the third derivative of these curves. 

Given the increased complexity of the phase behavior 
(separation of cloud and shadow curves; critical point not 
at the maximum of these; finite width coexistence region 
on the p—T plot collapsing to a line in the monodisperse 
limit), we find this somewhat remarkable. We emphasize 
that every term in Eq. (31) must be retained to achieve 
this degree of consistency - and in particular to reproduce 
the correct curvature of the cloud; to get the curvature 
of the shadow also correct, one needs as well the higher 
terms in Eq. (45). The latter were included in the figure 
although leaving them out gives a result that, for this 
example, is only slightly less accurate. Despite all the 
extra complexity, the region of the phase diagram within 
which the Landau expansion could tolerably be used to 
replace the full calculation is similar to the monodisperse 
case. 



which, unlike our work, it is not limited in principle). 
Finally, we note that although our Landau expansion 

was derived within the specific context of Eq. (9), the 
generic features of the phase diagram that it reproduces 
are seen rather generally in experiments and other theo- 
retical models. We have not attempted explicit calcula- 
tion for the case where more than one moment is present 
in the excess free energy, but since the topology of the 
phase diagram (near the fluid-fluid critical point) is al- 
ready captured by Eq. (9), it is possible that any such 
refinements will not alter the basic structure of these re- 
sults. A firm conclusion on this point is, however, left to 
future work. 

Acknowledgements We are grateful to Yuri Martfnez- 
Raton, Jose Cuesta, Patrick Warren and Peter Sollich for 
discussions. Work funded by EPSRC Grant GR/M29696. 
C.R. also acknowledges support by the Ministerio dc 
Ciencia y Tecnologia (Spain) under grant BFM2000- 
0004. 



VI. CONCLUSIONS 

Within the one-moment form of the excess free en- 
ergy represented by Eq. (9), we have derived a gener- 
alized Landau expansion, Eqs. (31,45), from which the 
characteristic features of polydisperse binary fluid coexis- 
tence can all be shown to follow. This Landau expansion 
was compared with direct calculation from the underly- 
ing (mean-field) free energy in a specific example (Fig- 
ures 1 and 2) with excellent agreement near the critical 
point and qualitative agreement away from it. Note that 
truncation of the Landau expansion at the quartic order 
(Eq.31) may prove adequate for qualitative purposes if 
the curvature of the shadow is not important. 

Our expansion is correct around the critical point for 
arbitrary polydispersity, be it broad or narrow. The 
monodisperse limit is recovered by decreasing the poly- 
dispersity. In a pressure-temperature representation this 
limit is highly singular: for example the slope of the 
perimeter of the coexistence region at the critical point 
will converge, as the monodisperse limit is taken, to a 
negative value — even though the coexistence line (onto 
which the region itself collapses in this limit) has posi- 
tive slope everywhere! This behavior is easily understood 
from Figure 2, and the Landau expansion seems to have 
no difficulty dealing with it. 

The perturbation series in density that underlies the 
Landau expansion contrasts with, and is complementary 
to, a previous perturbative approach that was developed 
for the limit of narrow polydisperse distributions [12]. 
The latter approach always breaks down close enough 
to the critical point (even within mean field theory, to 
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Pi/ Pi 



FIG. 1. Cloud (continuous) and shadow (dashed) curves of 
the polymeric system described in the text. The monodispcrsc 
limit, AN ^0, for which cloud and shadow curves merge into 
a single curve (the binodal), is also plotted as a dotted line. 
The thin lines represent the exact numerical solution of the 
model, while the thick lines represent the Landau expansion, 
which shows exact first and second order derivatives at the 
critical point. In the monodisperse case, the critical point 
rests at the top of the curve. 
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FIG. 2. p — T phase diagram of the polymeric system de- 
scribed in the text. The thin line represents the exact nu- 
merical solution of the model, while the thick lines represent 
the Landau expansion. This curve shows exact first and sec- 
ond order derivatives at the critical point. The monodisperse 
limit is also shown (dotted line). Inset: Magnified view of the 
critical point. 



Quantity 


Landau Coefficients 


Moments 


acL 


Ai,A2 


ni,n2,n3 


asH 


Ai,A2,Bi,B2,Co 


ni,n2,n3,n4 


PCL 


Ai-4, Bi,B2,Co 


ni,n2,n3,n4 


PSH 


Ai-4, Bi-4, Co,Ci,C2,Do 


ni,n2,n3,n4,n5 



TABLE I. Landau coefficients and moments needed to 
calculate slope and curvature parameters of the cloud and 
shadow. 
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